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ABSTRACT. The large genus Begonia began to diverge in Africa during the Oligocene. The 
current hotspot of diversity for the genus in China and Southeast Asia must therefore be the 
result of an eastward dispersal or migration across the Asian continent. To investigate the 
role of the Himalayas as a mesic corridor facilitating this migration, we constructed a time- 
calibrated molecular phylogeny using ITS sequence data. Himalayan species of Begonia were 
found to fall into two groups. The first is an unresolved grade of tuberous, deciduous species 
of unknown geographic origin, with evidence of endemic radiations in the Himalayan region 
beginning c. 7.4 Ma, coinciding with the onset of the Asian monsoon. The second is a group 
of evergreen rhizomatous species with a probable origin in China, which immigrated to the 
Himalayan region c. 5.1 Ma, coinciding with an intensification of the monsoon. The hypothesis 
of the Himalayas being a mesic migration route during the colonisation of Asia is not refuted, 
but further data is needed. 
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Introduction 

The large genus Begonia has around 750 species in Asia, with the bulk occurring 
in Southeast Asia, and the Malesian region having c. 440 of these (Hughes 2008). 
Reconstructions of the phylogenetic history of Begonia show an early divergent 
African clade, with Asian and American species nested within an African grade (Plana 
et al. 2004, Forrest et al. 2005, Goodall-Copestake et al. 2010) and Asian (including 
Socotran) species supported as monophyletic (Goodall-Copestake et al. 2010, Thomas 
2010). Given the initial diversification of Begonia in Africa in the Oligocene or late 
Eocene (Goodall-Copestake et al. 2009) and the monophyly of Asian species, there 
must have been a migration or dispersal of an ancestor eastwards from Africa which 
has speciated and led to the various diversity hotspots across Asia. 
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Molecular divergence age estimates indicate that the origins of Asian 
(including Socotran) Begonia date to 18-15 Ma (Thomas 2010) coinciding with the 
mid-Miocene climatic optimum, a warm phase which led to the expansion of tropical 
vegetation in Asia as far north as southern Japan and as far east as the northwest of the 
Indian subcontinent (Zachos et al. 2001, Morley 2007). However it would seem that a 
straightforward migration of Begonia across the Arabian peninsula to Asia is unlikely, 
given the dry conditions thought to have existed across much of southwest Asia and 
the Arabian peninsula at the time (Morley 2007). Hence a long distance dispersal event 
does not seem an unlikely hypothesis for the entry of an African ancestor into Asia, 
although this was possibly facilitated by the greater expanse of tropical climate during 
the mid-Miocene. Begonia from Sri Lanka, the Western Ghats and Socotra are at the 
western limit of the Asian Begonia clade (Fig. 1), and have been found to be early 
divergent lineages within a phytogeny of Asian species (Thomas 2010). 

The Himalayas began to form during the initial collision of the Asian and 
Indian plates during the early Eocene c. 35 Ma (Ali & Aitchison 2008) and by the time 
of the arrival of Begonia in Asia c. 15 Ma, the High Himalayas were present (Amano 
& Taira 1992). Hence, depending on climate, there was potentially a significant area 
of montane habitat at the north of the Indian continent which could have provided 
a mesic link as a migration route eastwards towards the current Begonia diversity 
hotspots of lndo-China and Malesia (Fig. 1). 

To investigate the role of the Himalayas as an easterly migration route during 
the colonisation of Asia, Himalayan species of Begonia were placed in a phytogeny 
to allow their interpretation in a geographical and temporal context. This paper builds 
on previous studies by considerably increasing the sampling of Himalayan species 
and coding the region as a distinct geographic unit in a biogeographic analysis. In 
particular, this paper addresses the following questions: 

Is there evidence for the Himalayas behaving as a link in an eastward migration of 
Begonial 

—How has paleo-climatic change and the Himalayan orogeny influenced the evolution 
of Nepalese Begonial 


Materials and methods 

The molecular phylogenetic analysis was based on nuclearribosomal internal transcribed 
spacer (ITS) sequences, obtained from Genbank and also de novo for this study 
(Appendix A) following the methods in Thomas (2010). Phylogenetic reconstruction, 
divergence times and ancestral area reconstruction were performed simultaneously 
using Bayesian inference as implemented in BEAST vl.6.1 (Drummond & Rambaut 
2007), using 4 x 10 7 generations and sampling every 1000th generation. The sequence 
data was divided into three partitions, namely the two internal transcribed spacers 
(ITS1 and ITS2) and the 5.8s ribosomal gene. Models of sequence evolution for each 
nucleotide sequence partition were determined using jModelTest (Posada 2008). The 
divergence time for the split between African and Asian species was set as 17.8 Ma 
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Fig. 1. The distribution of Begonia in the Old World, highlighting the disjunction between the 
African and Asian floras, and the westward linear distribution along the Himalayas. Data from 
the GBIF data portal and Hughes & Pullan (2007). 


(Thomas 2010). implemented as a normally distributed prior with a standard deviation 
of 3.4. Bayesian ancestral area reconstructions were performed using the continuous¬ 
time Markov chain (CTMC) model for discretized diffusion specified by Lemey et al. 
(2009). considering diffusion amongst 10 geographic areas defined as Africa. China. 
Himalaya. India, Indo-China. Philippines. Socotra. Sulawesi. Sunda Shelf and Taiwan. 
Maximum clade credibility consensus trees were made using TreeAnnotator v 1.6.1 
(Drummond & Rambaut 2007) with a bum-in of 25%. and visualised using FigTree 
vl.3.1 (Rambaut 2007). 

The Bioclim parameter precipitation seasonality (Busby 1991) was extracted 
from the WorldClim database (Hijmans et al. 2005a) using DIVA-GIS (Hijmans et al. 
2005b) and georeferenced herbarium specimen data from Hughes & Pullan (2007). It 
is calculated as the standard deviation of weekly precipitation estimates expressed as 
a percentage of the annual mean of those estimates. Scores for species with multiple 
localities were averaged. 

Perennation organ type was scored from literature records and herbarium 
specimens available from Hughes & Pullan (2007). 
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Results 


Sequence characterisation 

The aligned sequences gave a dataset of 873 characters in length (1TS1 330 bases, 
5.8s 148 bases, ITS2 395 bases). The bases at positions 492^499 and 505-525 were 
excluded from the analysis due to difficulty in reliably assigning positional homology. 
The remaining 845 bases were included in the analysis, of which 561 were variable and 
433 (51%) were potentially informative (ITS 1,237 variable and 190 (58%) potentially 
informative; 5.8s, 31 variable and 11 (7%) potentially informative; ITS2, 294 variable 
and 232 (63%) potentially informative). All three regions were analysed under a 
GTR+G+I model (general time reversible model, gamma distributed rate variation, 
plus a proportion of invariable sites). 

Phylogenetic analysis—geography and dating 

The Asian + Socotra ingroup was supported as monophyletic with a posterior 
probability of 1 (Fig. 2). Relationships at the base of the tree were generally poorly 
supported, and hence inferences regarding timing and geography of cladogenesis on 
this part of the tree are not possible with confidence. 

Clade A contains species from China, Indo China and Malesia, with several 
well supported sub-clades. One of these contains the Thai species B. alicida and 
B. smithiae , together with three species of Begonia sect. Pan'ibegonia. Two other 
supported subclades are highlighted (Fig. 2) concerning the large Begonia sect. 
Petermannia with a probable origin on the Sunda Shelf in the late Miocene followed 
by later diversification in Wallacea. 

Himalayan species are present only in Clade B (Fig. 2), in both the 'Diploclinium 
grade’ and the clade consisting of members of Begonia sect. Platycentrum and 
Sphenanthera , marked PLA-SPH in Fig. 2. The 'Diploclinium grade’ has a largely 
unsupported topology at the base, but began diversifying sometime between 14.7 
(±7.6) and 12.3 (±6.7) Ma in the mid-Miocene. Species from the Himalayan region 
in this grade are intermixed with species from China, Indo-China and Taiwan. Due 
to the unresolved backbone, the area of origin for this tuberous, seasonally adapted 
grade cannot be reconstructed. Supported subclades within the grade have areas 
of origin reconstructed as the Himalaya, Indo-China or China. One consists of the 
Taiwan endemic B. ravenii sister to two accessions of the Himalayan B. dioica ; the 
geographic origin of this clade is equivocal, with the Himalayas having the highest 
posterior probability, though only of 0.37; it dates from 7.4 (±5.7) Ma. Another clade 
consists of four species from Thailand, probably representing speciation originating 
in that country (PP 0.84), dating from 8.1 (±5.2) Ma. A clade of species from Nepal, 
with B. picta appearing as sister to the remaining species, likely has a Himalayan 
origin (PP 0.90) dating 7.4 (±4.6) Ma. A clade with Begonia putt'd (Thailand) and B. 
labordei (China) along with an unidentified species from China diverged 6.4 (±4.8) 
Ma in China (PP 0.82). The most highly nested subclade in the Diploclinium grade, 
sister to the PLA-SPH clade, contains three accessions of Begonia rubella from Nepal 
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Fig. 2. Maximum clade credibility chronogram derived from an analysis of nuclear ribosomal 
ITS sequences using BEAST v 1.6.1 (Drummond & Rambaut 2007). Asterisks denote 
clades with a posterior probability of less than 0.8; pie charts show Bayesian ancestral 
area reconstructions on supported nodes; lower scale denotes time in millions of years; the 
monochrome and coloured boxes at the branch tips show the perenniation organ type and 
precipitation seasonality respectively; missing data coded as white. Three-letter codes indicate 
the sections of Begonia (ALI, sect. Alicida; AUG, sect. Augustia; COE, sect. Coelocentnnn; 
DIR sect. Diploclinium; HAA, sect. Haagea ; IGN, ignota (section unknown); MON, sect. 
Monopteron ; PAR, sect. Panibegonia; PET, sect. Petermannia; PLA, sect. Platycentmm ; REI, 
sect. Reichenheimea; RID, sect. Ridleyella; SPH, sect. Sphenanthera). 
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and one of B. tribenensis , also from Nepal. This subclade probably has an origin in the 
Himalayas (PP 0.95) dated 6.4 (±4.3) Ma. 

The origin of the PLA-SPH clade dates from 6.5 (±4) Ma and hence is 
considerably younger than the 'Diploclinium grade’. Himalayan species in the PLA- 
SPH clade are intermixed with species from China, Indo-China, the Sunda shelf, 
Taiwan, the Philippines and Sulawesi. The geographic area of origin for this clade is 
most probably China (PP 0.73). Three supported subclades within the PLA-SPH clade 
also have an area of origin optimised as China, with posterior probabilities of 0.65, 
0.70 and 0.94. Only one subclade has a probable area of origin in the Himalaya (PP 
0.99), dating from 3.0 (±2.3) Ma and with one dispersal into Indo-China. The overall 
picture for the PLA-SPH clade is of a fairly rapid diversification with its origins in 
China following loss of the tuberous habit, with dispersal to the Himalayan region and 
throughout Southeast Asia followed by localised speciation. 

Phylogenetic analysis—perennation organs , precipitation and seasonality 
All of the species in the "Diploclinium grade’ are tuberous and the vast majority are 
found in climates with a marked seasonality in rainfall (Fig. 2). They tend to die back 
and lose their leaves completely in the dry season. Species within the PLA-SPH clade 
either have a rhizome or no specialised perennation organs. The rainfall seasonality 
of the distributions of most of the species is much less marked than for those in the 
‘Diploclinium grade’, and they do not lose their leaves during the dry season. Only 
one clade lacks perennation organs entirely, consisting of Begonia sect. Petermannia 
species distributed in the largely everwet Malesian region. 


Discussion 

Although support for the backbone of the tree was weak, the taxon composition 
of the two major clades, A and B, matches that of strongly supported clades found 
using chloroplast sequence data by Thomas (2010), with much greater sampling of 
Himalayan species in this study. In addition, the overall topology of clade B matches 
the results of Thomas (2010) with respect to a clade of Begonia sects. Platycentrum 
and Sphenanthera (PLA-SPH in Fig. 2) being nested within a grade consisting mainly 
of Begonia sect. Diploclinium species from the continent (labelled ‘Diploclinium 
grade’ in Fig. 2); also congruent is the Chinese species B. grandis being sister to the 
remaining taxa. The basal relationships in Clade A are unsupported, but the supported 
subclades within Clade A are similar to those found by Thomas (2010), and generally 
correspond to species of a single section or geographic region. This congruence gives 
some confidence in the tree topology, despite the disappointing levels of support. 

The Himalayas as a link in an eastward Begonia migration 

If the Himalayas acted as a migration corridor in the early evolution of Asian Begonia , 
we would expect to see other species in Asia nested within a Himalayan grade in the 
phylogeny. Due to the unsupported nodes at the base of the clade containing all the 


Origin and divergence of Asian Begonia 


283 


Himalayan species, evidence for this is equivocal. Of the two major sections in the 
Himalayan Begonia flora {Begonia sects. Diploclinium and Platycentmm ). Begonia 
sect. Diploclinium has the oldest lineages which date to between 14.7 (-7.6) and 12.3 
(±6.7) Ma and accounts for nearly 60% of the species in Nepal. Some subclades in 
this group do show a likely Himalayan origin, but the supported nodes are dated in the 
late Miocene at c. 7.5 Ma and so post-date the early divergence of the section. Further 
work on other gene regions is needed to resolve the relationships in the ‘Diploclinium 
grade'. The hypothesis of basally branching Himalayan lineages in Begonia sect. 
Diploclinium is not refuted by our results. However, the fact that the Chinese species 
B. grandis is possibly sister to the rest of the 'Diploclinium grade’ is tantalising and 
raises the possibility of China being an ancestral area for this group. 

The other major section. Begonia sect. Platycentnim. accounts for 35% of 
Nepalese species. The species are all included in the PLA-SPH clade with an origin 
in the late Miocene-earlv Pliocene. There is no evidence for the Himalayan species 
in Begonia sect. Platycentrum being relicts from an eastward migration early in the 
evolution of Asian Begonia , in fact the emerging picture is the exact opposite. This 
clade has a probable Chinese origin, and represents an entry of the genus into the 
Himalayan region separate from those in the ‘Diploclinium grade'. The Himalayan 
species are largely highly nested within this clade and there is evidence for more 
than one migration into the Himalayan region, including a very recent one in the late 
Pliocene-early Pleistocene (bottom of Fig 2: clade from B. hatacoa to B. nmvakotensis). 

Other major clades containing the bulk of species diversity in Southeast 
Asia ( Begonia sect. Petennannia from Malesia: Begonia sect. Diploclinium from the 
Philippines: Begonia sect. Reichenheimia from the Sunda shelf) remain unsupported 
in Clade A. Our results do not show any affinity of these clades with Himalayan taxa. 
and the origin of these large radiations remains enigmatic. 

Paleo-climatic change , Himalayan orogeny and evolution of Xepalese Begonia 
It is possible to highlight two main events in the evolution of the Himalayan Begonia 
flora—the mid- to late Miocene diversification of the tuberous, seasonally adapted 
‘Diploclinium grade' and the late Miocene-Pliocene immigration and diversification 
of the rhizomatous and evergreen PLA-SPH clade. 

The early diversification of tuberous clades with a likely origin in the 
Himalayas can be dated to between 6.4 (=4.8) and 7.4 (=4.6) Ma. This coincides 
with the development of the Asian monsoon 7.4 Ma (Copeland 1997), due to the 
Tibetan plateau having reached sufficient altitudes to affect a major re-organisation 
of atmospheric circulation over the Asian continent (Zheng et al. 2004). The tuberous 
species would be well adapted to this seasonal monsoon climate, with concentrated 
periods of intense rainfall interspersed with significant dry seasons. These species are 
deciduous, and their flowering periods are strongly constrained to the three months 
after the onset of the monsoon (Rajbhandary et al. 2010). The subsequently formed 
dry dehiscent capsular fruits are then able to disperse their seeds either passively or 
through wind assistance during the following dry season. 
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The immigration of members of the PLA-SPH elade into the Himalayan region 
began 5.1 (±3.2) Ma, with an endemic radiation starting 3.0 (±2.1) Ma and continuing 
throughout the Pleistocene. The immigration coincides with a further intensification of 
the Asian monsoon around 5 to 3.6 Ma, possibly linked to a further surge in the uplift 
of the Himalayas and the Tibetan plateau (Zheng et al. 2004, Zhisheng et al. 2001) 
and changes in ocean currents in the lndo-Pacific region (Srinivasan & Sinha 2000). 
The evolution of an evergreen, rhizomatous habit in this clade suggests adaptations to 
wetter conditions, but the details of how the Pliocene monsoon intensification affected 
total and seasonal rainfall in the Chinese (presumably SW China) region where the 
PLA-SPH clade originated are unknown. The Himalayan species in this clade currently 
occur in areas with similar rainfall seasonality 7 to the tuberous species, but occupy 
different micro habitats where water is more constantly available; their flowering is 
not strongly constrained by the monsoon (Rajbhandary et al. 2010). They have fruits 
which are adapted to rain-splash seed dispersal, and hence depend on either rainfall 
or drips and splashes from waterfalls to disperse. Our results are congruent with those 
of Tebbitt et al. (2006), who suggested that members of Begonia sect. Platycentrum , 
which have fruit morphologies indicative of rain dispersal, evolved from wind- 
dispersed Asian taxa following the colonisation of wetter habitats. 

Of the two groups of Begonia in the Himalayas, seasonally deciduous 
and evergreen, only the former could possibly be relicts of an eastward migration. 
However, due to the unresolved relationships of clades within this group, their 
origins remain unknown. The evergreen species represent a re-entry to the Himalayan 
region, most likely from China, and have speciated in response to further tectonic 
uplift and Pliocene-Pleistocene climatic cycles and changes in the monsoon intensity. 
These factors are also likely to have strongly influenced the recent diversification of 
other elements of the Himalayan flora. Further phylogenetic studies of Himalayan 
plants will reveal the relative contribution of relict clades and recent speciation to the 
considerable plant diversity in the region. 


ACKNOWLEDGEMENTS. We are grateful to the M.L. MacIntyre Trust and the University 
Grants Commission (UGC) Nepal for funding, to Ching-1 Peng for providing DNA sequence 
data, and to Michelle Hollingsworth, Alex Clarke and Ruth McGregor for assistance in the lab. 
This research was supported by the Scottish Government’s Rural and Environment Research 
and Analysis Directorate. 


References 

Ali, J.R. & Aitchison, J.C. (2008) Gondwana to Asia: Plate tectonics, paleogeography 
and the biological connectivity of the Indian sub-continent from the Middle 
Jurassic through latest Eocene (166-35 Ma). Earth Sci. Rev 88: 145-166. 
Amano, K. & Taira, A. (1992) Two-phase uplift of Higher Himalayas since 17 Ma. 
Geology> 20: 391-394. 


Origin and divergence of Asian Begonia 


285 


Busby, J.R. (1991) BIOCLIM - A Bioclimatic Analysis and Prediction System. In: 
Margules, C.R. & Austin, M.P. (eds) Nature Conservation: Cost Effective 
Biological Surveys and Data Analysis, pp. 64-68. Canberra: CSIRO. 

Copeland, P. (1997) The when and where of the growth of the Himalaya and the 
Tibetan Plateau. In: Ruddiman, W.F. (ed) Tectonic Uplift and Climate Change. 
New York: Plenum Press. 

Drummond, A.J. & Rambaut, A. (2007) BEAST: Bayesian evolutionary analysis by 
sampling trees. BMC Evol. Biol. 7: 214. 

Forrest, L.L., Hughes, M. & Hollingsworth, P.M. (2005) A phylogeny of Begonia 
using nuclear ribosomal sequence data and morphological characters. Svst. Bot. 
30: 671-682. 

Goodall-Copestake, W.P., Harris, D.J. & Hollingsworth, P.M. (2009) The origin of 
a mega-diverse genus: dating Begonia (Begoniaceae) using alternative datasets, 
calibrations and relaxed clock methods. Bot. J. Linn. Soc. 159: 363-380. 

Goodall-Copestake, W.P., Perez-Espona, S., Harris, D. & Hollingsworth. P.M. (2010) 
The early evolution of the mega-diverse genus Begonia (Begoniaceae) inferred 
from organelle DNA phylogenies. Biol. J. Linn. Soc. 101: 243-250. 

Hijmans, R.J., Cameron, S.E., Parra, J.L., Jones, P.G. & Jarvis, A. (2005a) Very high 
resolution interpolated climate surfaces for global land areas. Intern. J. Cliniatol. 
25: 1965-1978. 

Hijmans R.J., Guarino, L., Mathur, P. & Jarvis, A. (2005b) DIVA-GIS Version 5.0. 
http://www.diva-gis.org/ 

Hughes, M. (2008) An Annotated Checklist of Southeast Asian Begonia. Edinburgh: 
Royal Botanic Garden Edinburgh. 

Hughes, M. & Pullan, M. (2007) Southeast Asian Begonia Database. Electronic 
publication accessible via www.rbge.org.uk. 

Lemey, P., Rambaut, A., Drummond, A.J., & Suchard. M.A. (2009) Bayesian 
phylogeography finds its roots. PLoS Computat. Biol. 5. 

Morley, R.J. (2007) Cretaceous and Tertiary climate change and the past distribution of 
megathermal rainforests. In: Flenley, J.R., & Bush, M.B. (eds) Tropical Rainforest 
Responses to Climatic Changes , pp. 1-31. Berlin: Springer. 

Plana, V., Gascoine, A., Forrest, L.L., Harris, D. & Pennington, R.T. (2004) Pleistocene 
and pre-Pleistocene Begonia speciation in Africa. Molec. Phylogen. Evol. 31: 
449-461. 

Posada, D. (2008) jModelTest: Phylogenetic model averaging. Molec. Biol. Evol. 25: 
1253-1256. 

Rajbhandary, S., Shrestha, K.K. & Hughes, M. (2010) Distribution patterns of Begonia 
species in the Nepal Himalaya. Botanica Orientals 7: 73-78. 

Rambaut, A. (2007) FigTrce, software package downloaded from http://tree.bio.ed.ac. 
uk/software/figtree/ 

Srinivasan, M.S. & Sinha, D.K. (2000) Ocean circulation in the tropical Indo-Pacific 
during early Pliocene (5.6-4.2 Ma): Paleobiogeographic and isotopic evidence. 
Proc. Indian Acad. Sci., Earth Planet. Sci. 109: 315-328. 


286 


Card Bull. Singapore 63(1 & 2) 2011 


Tebbitt, M.C., Forrest, L.L., Santoriello, A., Clement, W. & Swensen, S.M. (2006) 
Phylogenetic relationships of Asian Begonia , with an emphasis on the evolution 
of rain-ballist and animal dispersal mechanisms in sections Platyeentrum , 
Sphenanthera and Leprosae. Syst. Bat. 31:327-336. 

Thomas, D.C. (2010) Phylogenetics and Historical Biogeography of Southeast Asian 
Begonia (Begoniaceae). U.K.: University of Glasgow. 

Zachos, J., Pagani, M., Sloan, L., Thomas, E. & Billups, K. (2001) Trends, rhythms 
and aberrations in global climate 65 Ma to present. Science 292: 686-693. 

Zheng, H., Powell, C. McA., Rea, D.K., Wang, J. & Wang, P. (2004) Late Miocene and 
mid-Pliocene enhancement of the East Asian monsoon as viewed from the land 
and sea. Global Planet Change 41: 147-255. 

Zhisheng, A., Kutzbach, J.E., Prell, W.L. & Porter, S.C. (2001) Evolution of Asian 
monsoons and phased uplift of the Himalaya-Tibetan plateau since late Miocene 
times. Nature 411: 62-66. 


Appendix A. List of the Genbank accession numbers of the DNA sequences used in the 
phylogenetic analysis. 

B. aceroides, HQ729023; B. acetosella, AY048976; B. aequata, AF48514; B. alicida , HQ729022; 
B. alpitia, AY753717; B. alveolate, AY048977: B. ampliioxus, AF485150; B. annulate, IIQ729060; 
B. annulate2, HQ729059; B. aptera, AF48510; B. atricha , HQ729047; B. balansana, AF485091; 
B. barbate, AF485095; B. baviensis, AF485087; B. brevipes, HQ729048; B. biyophila, HQ729030; 
B. bnrbidgei, HQ729049; B. eathavana, AF280106; B. ceratocarpa, AY048978; B. cliiastnogy'tia, 
HQ729050; B. cliloroneura, AF485134 ; B. chiorostricta, AF485153; B. circnrnlobata, AY753721; 
B. cleopatrae, AF485133: B. crispipila, HQ729051; B. demissa, HQ729026;/?. dioical, IIQ729039; 
B. choice 2, HQ729038; B. dipetala, AF485142; B . discrete, HQ729024; B. dregei, AY429336; 
B. flagellar is, HQ729031; B.fonnosana, AF485119; B. geranioides, AF469120; B. goegoensis, 
AF485138; B. grandis, AF485089; B. grandisJ, AF485088; B. halconensis, AF485106; B. handelii, 
AY048982; B. hatacoal, AF485111; B. hemsleyana, AF485099; B.hernandioides, AF485135; B. 
livatae, AJ287262; B. incise, AF485148; B. Josephii2, HQ729037; B. kingiana, AF485139; B. 
koordersii, HQ729052; B. labotxlei, AF485122; B. laruei, HQ729058; B. leprosa, AY753722; B. 
leptoptera, HQ729036; B. leptoptera2, HQ729033; B. feptoptera3, HQ729034; B. leptoptera4, 
HQ729035; B. longicarpa, AF485109; B. longifolia, AF485105; B. macintyreana, HQ729054; 
B. tnalabarica, AF468141; B. malacbosticta, AF485156; B. masoniana, AF485123; B. morsei, 
AF485130; B. multangnla , AY753724; B. mnricata, AY753725; B. negrosensis, 11Q729055; 
B. nepalensis, AY753726; B. nmvakotensis, 1IQ729061; B. ovatifolia, HQ729032; B. palmate, 
AF485113.1; B. palmate, AF485115.1; B. panclitliarensis, HQ729062; B. parvula, GUI76066; B. 
pasamanensis, HQ729070; B. pavonine, AY753727; B. pitta, HQ729042; B. picta2, HQ729041; 
B. psendolateralis, HQ729053; B. putii, HQ729025; B. rabillii, HQ729027; B. rajah, AF485136; 
B. ramosii, HQ729057; B. ravenii, HQ729040; B. rex, AY753728; B. robnsta, AY753729; B. 
roxburghii, AF485092; B. rubella, AF485112; B. rnbella2, HQ729043; B. mbella5, HQ729044; 
B. rnboides, AY048987; B. samliaensis, AF469122; B. scottii, HQ729063; B. setifolia, AY048990; 
B. sikkimensis, HQ729064; B. silletensis, AF485094; B. sinwaensis, HQ729029; B. smitliiae, 
HQ729065; B. socotrana, AF469121; B. snborbiculata, HQ729069; B . taiwaniana, AY753731; 
B. taligeral, HQ729066; B. taligera2, HQ729067; B. ten uifolia, HQ873478; B. teysmanniana, 
HQ729068; B. tribenensis, HQ729045; B. variabilis, AY753732; B. variabilis, HQ729046; B. 
varipeltata, HQ729056; B. versicolor, AF485090; B. wenslianensis, AY048974; B. xantliina, 
AY753733; Begonia sp. (China), GUI76063; Begonia sp. (Nepal). HQ729028. 



